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Abstract: We calculate the cosmic ray muon flux at ground level using directly the primary cosmic ray 
spectrum and composition measured in the ATIC-2 balloon experiment. In order to extend the calculations 
to more high energies, up to 100 TeV, we use the data of the GAMMA experiment as well as Zatsepin 
and Sokolskaya model. This model supported by the ATIC-2 data comprises contributions to the cosmic 
ray flux of three classes of astrophysical sources - the shocks from exploding stars, nova and supernova 
of different types. The muon flux computation is based on the method for solution of atmospheric hadron 
cascade equations in which rising total inelastic cross-sections of hadron-nuclear interactions as well as 
non-power-law character of the primary cosmic ray spectrum are taken into account. The calculated muon 
spectrum agrees well with measurements of L3+Cosmic and BESS-TeV, CAPRICE, Frejus, MACRO, 
LVD as well as other experiments. 
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Introduction 

The muons produced through the cosmic ray in- 
teractions with the Earth atmosphere provide the 
tool for indirect study of the primary cosmic ray 
(PCR) spectra. May comparison of the predicted 
and measured atmospheric muon (AM) flux serve 
as reliability trial for PCR data? The answer de- 
pends on the relationship between size of the PCR 
uncertainties and that of AM flux. To attempt an- 
swering the question we calculate the cosmic ray 
muon flux at the ground level using directly the 
data on PCR spectrum and composition measured 
in the ATIC-2 experiment [1]. 

In order to compare the predictions with the high- 
energy measurements of the AM flux we extend 
the calculations to more high energies, up to 100 
TeV, using also the PCR spectrum data of the 
GAMMA experiment [2]. The PCR model by Zat- 
sepin and Sokolskaya [3] supported by the ATIC- 
2 data was applied as the nice instrument to ex- 
trapolate median energy data to high energy one. 
This model comprises contributions to the cosmic 
ray flux of three classes of astrophysical sources 
like supernova and nova blast waves (shocks). The 
muon flux calculation is based on the method 



to solve the atmospheric hadron cascade equa- 
tions [4, 5] in which we take into account rising to- 
tal inelastic cross-sections of hadron-nuclear inter- 
actions as well as non-power-law character of the 
primary cosmic ray spectrum. A high convergence 
of the method provides an operative way to calcu- 
late the secondary cosmic ray fluxes and allows to 
test "on-the-fly" the primary spectrum models. 

Hadron fluxes were computed with slightly revised 
Kimel and Mokhov parametrization (see [6, 4]) 
for nucleon and meson production cross sections 
which are close to the SIBYLL mini-jet model [7]. 

Primary cosmic ray spectra 

The balloon borne experiment ATIC (Advanced 
Thin Ionization Calorimeter) [1], designed for 
measurements of cosmic rays energy spectra with 
individual charge resolution from protons to iron, 
enabled to obtain PCR spectra in the wide energy 
interval 50 GeV - 200 TeV with high statistical as- 
surance. The differential spectra of protons and he- 
lium nuclei obtained in the ATIC-2 experiment are 
shown in figure 1 along with a bulk of data from 
balloon, satellite and ground based experiments - 
BESS [8], AMS [9], IMAX [10], CAPRICE [11], 
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Figure 1 : Primary proton and helium spectra, combining balloon, satellite and ground-based measurements. 
The solid curve presents Zatsepin and Sokolskaya model [3]. Eq is the kinetic energy of the particle. 



MASS [12], RICH [13], MUBEE [14], RUN- 
JOB [15], JACEE [16], SOKOL [17], KASCADE 
SH [18], GAMMA [2], HEGRA [19], TIBET 
HD [20], ICHIMURA [21]. 

Proton and helium spectra measured in the ATIC-2 
experiment have different slopes and differ from 
a simple power law. The ATIC-2 data are in 
agreement with the data of magnetic spectrome- 
ters (BESS, AMS, IMAX, CAPRICE, MASS) be- 
low 100 GeV. In the energy region 1 < E < 
10 TeV the ATIC-2 data are consistent with the 
SOKOL measurements and with those of atmo- 
spheric Cherenkov light detector HEGRA. At en- 
ergies above ^10 TeV the spectra become steeper, 
and follow the data of emulsion chamber exper- 
iments MUBEE and JACEE, though the agree- 
ment is not so clear. The solid curves in figure 1 
are to present the model suggested by Zatsepin 
and Sokolskaya (ZS) [3] that fits well the ATIC- 
2 experimental data and describe PCR spectra in 
the energy range 10-10 7 GeV. In order to extend 
our calculation to higher energies, the PCR spec- 
tra measured in the GAMMA [2] experiment was 
used. The energy spectra and elemental composi- 
tion, obtained in the GAMMA experiment cover 
the 10 3 -10 5 TeV range (shaded areas) and agree 
with the corresponding extrapolations of known 



balloon and satellite data at the E > 10 3 TeV. 
In the present calculation, a version of the spec- 
tra, reconstructed in the framework of 1 . 2D com- 
bined analysis with the SIBYLL interaction model 
(see [2] for details), was utilized. 

Conventional atmospheric muons 

Apart from evident sources of AM, ir M 2 and 
decays, we take into consideration three -particle 
semileptonic decays, K^ 3 , K® 3 . Also we take into 
account small fraction of the muon flux originated 
from decay chains K — > tt — > fi (Kg — > 7r + +7r~, 

K± -> TT* + 7T°, Kl -> 7T ± + + U t {v t ), 

t = e, fi). We do not consider here a conjec- 
tural prompt muon component of the flux (see e.g. 
[22, 23]). In figures 2, 3 presented are results of the 
calculation of the surface muon flux along with the 
data of muon experiments that comprise the direct 
measurements of CAPRICE [24], BESS-TeV [25], 
L3+Cosmic [26], Cosmo-ALEPH (see Ref. [27]), 
L3 and MASS (taken from [22]) as well as the 
data (converted to the surface) of underground ex- 
periments MSU [28], MACRO [29], LVD [30], 
Frejus [31], Baksan [32], Artyomovsk [33]. The 
light shaded areas in figure 2 and figure 3 (the left 
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Figure 2: Energy spectrum of muons at ground level near vertical. The dashed-line curves and the shaded 
area present this work calculation with the ATIC-2 primary cosmic-ray spectrum. 
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Figure 3: High-energy plot of the ground level muon spectrum. The dashed-line curves and shaded areas 
present this work calculations with the ATIC-2 primary spectrum (left) and GAMMA one (right). The solid 
curve presents the calculation with Zatsepin and Sokolskaya model. 
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corner) show the muon spectrum calculated with 
the ATIC-2 primary spectra taking into consider- 
ation statistical errors (dashed curve corresponds 
to mean values). For the range 10 — 3000 GeV 
one sees fair accordance of the muon flux, calcu- 
lated with the ATIC-2 spectra, and the recent mea- 
surements but the Cosmo-ALEPH data. The high- 
energy part of the muon flux is shown in figure 3, 
where the dark shaded area (at the right) presents 
our calculation with the GAMMA primary spectra 
input and the solid curve presents the muon flux 
computed with ZS primary spectrum model which 
appears to be a reliable bridge from TeV range to 
PeV one. It should be noted that without consid- 
ering the prompt muon contribution above 10 TeV 
one can say about satisfactory agreement of cal- 
culated fluxes only with the data of MACRO and 
LVD measurements. 

In conclusion, it may be said that the high accuracy 
of the ATIC-2 data results in the muon flux calcula- 
tion uncertainty, comparable with rather high pre- 
cision of the last decade muon flux measurements. 
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